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\mJ \ Abstract 

*vj ■ We address the formation of correlations in a general nonequilibrium 

^Nl ' situation. As an example we calculate the formation of light composite 

particles in a heavy ion collision. In particular we study the formation 
of deuterons via three-body reactions in some detail. To calculate the 
^ ' relevant reaction rates we solve the Faddeev equation which is consis- 

tently modified to include the self energy corrections and Pauli blocking 

^Q ■ factors due to surrounding nucleons. We find that the time scales of the 

reaction as well as the number of emitted deuterons change as medium 
dependent rates are used in the calculation. 



OS 

OS ' 1 Introduction 

The description of the dynamics of an interacting many-body system is partic- 
ularly difficult when the quasiparticle approach reaches its limits. That may 
be the case when the residual interaction is strong enough to build up corre- 
lations. An example of such a system with correlations is nuclear matter. In 
heavy ion collisions, the nuclear matter is first excited and compressed and 
then decompressed. On a macroscopic scale, the nuclear matter is formed dur- 
ing a supernova collapse and becomes the material of which a neutron star is 
V^j . made. 

C^ ' Here, we address the formation of correlations in a general nonequilibrium 

situation. The simplest process in nuclear matter is the formation of deuterons; 
their number is, besides the numbers of other composite particles (tritons, 
helium-3's, a-particles), an important observable of heavy ion collisions. A 
microscopic approach to treat this complicated process uses the Boltzmann 
equation. In this context the equation has been numerically solved utilizing 
the Beth-Uehling-Uhlenbeck approach, see e.g. Ref. p. For nucleon (/jy) and 
deuteron (/<j) Wigner distribution functions, the Boltzmann equation reads 

dtfN + {U, In} = JCjv[fN,fd] (1 - In) - ^'[/jVj/d] In, 
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dtf d + {uj d } = }Cf[f N j d ](i + f d )-)CT t [fN,Mfd, (i) 

where U is a mean-field potential and {•, •} denotes the Poisson brackets. The 
collision integrals JC appearing in Eq. (|l]) couple the distribution functions and 
respect all the collisions (elastic and reactive) between the constituents. We 
focus on the reaction part only, e.g. the loss reaction for the deuteron is given 
by 
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(2) 



where / = 1 ± / for bosons/fermions. The transition matrix element is for 
the break-up operator Up. and given by the solution of the Alt-Grassberger- 
Sandhas (AGS) equationcl explained below. Several approximations have been 
used to describe the transition matrix Uq, e.g. Born approximation that may be 
justified in the context of weak interactions, or impulse approximation that is 
justified for higher scattering energies. One further strategy is to introduce the 
cross section in the above equations and use experimental values or meaningful 
extrapolations (assuming detailed balance for the back reaction). 

2 In-medium three-body reactions 

As nuclear matter provides a dense system we address the question to what 
extend the reaction depends on the embedding medium. To this end we have 
derived an AGS-type equation for the three-particle correlation embedded in 
an uncorrelated mediumB&El. This will be shortly sketched here. The equation 
is derived in the context of the cluster expansion or Dyson equation approachO. 
The expansion is consistent in the sense that the respective in-medium one- and 
two-body problems are solved and implemented in the three-body equation, 
i.e. nucleon self-energies e, Pauli blocking factors, and the Mott effect of the 
deuteron are all consistently included. The cm. momentum is treated using 
four-body kinematics. Note, that the cm. momentum enters only parametri- 
cally into the three-body equation. The three-body Green function evaluated 
in an uncorrelated medium is then given by 
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Figure 1: Neutron deuteron crosn sections. Figure 2: Break-up cross section for differ- 
Experimental data Schwarz et al. a ent densities at a temperature T = 10 MeV. 



where the corresponding single particle self-energy in Hartree-Fock approxi- 
mation is given by 
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£ ffF (l), E ffF (l) = ^[y 2 (12,12)--/ 2 (12,21)]/ 2 , (4) 



We assume the dominant two-body interaction Vi only. For..a recent and 
throughout investigation of nuclear three-body forces see Ref. H. A similar 
equation has been derived for the channel] Green function G 3 (z), for the 
two-body subchannels. The AGS transition operator U a p(z) is defined by 



G 3 (z) = 8 afj G { «\z) + G[ a \z)U a0 {z)G<£\z). 



(5) 



For / — > this operator leads to the isolated AGS-operator with the correct 
reduction formulas. The resulting in-medium AGS equation for U a p(z) is then 
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To simplify notation we have introduced A3 = /1/2/3 + /1/2/3 and for the two- 
body subchannel of the three-body system A 2 = /1/2 — /1/2 (cycl. perm.) as 

tThe three-body fragmentation channels are labeled by the spectator particle a £ {1,2,3} 
and the break-up channel by zero. 



well as the three-body resolvent R 3 ' [z) = (z — e\ — £2 — £3) 1 - The respective 
two-body t matrix in the subchannel 7 entering in the above equation is 

T<?\z) = V^ + V^R^(z)N^T^(z). (7) 

The optical theorem has been used to calculate the break-up, the total and 
the elastic cross sections. For the isolated three-body case the respective the- 
oretical results are shown in Fig. [j] and compared to the latest experimental 
data on neutron deuteron scattering. The calculation is done using a rank one 
Yamaguchi potential for the 1 So and coupled 3 Si — 3 Di channels. The effect of 
the nuclear medium is seen in Fig. 0. The specific effects on the cross section 
are indicated in the figure. Note that the threshold shifts to smaller energies as 
it should, since the deuteron binding energy becomes smaller with increasing 
density. As the binding energy of the two-body subsystem approaches zero 
one might expect an infinite number of states in the corresponding three-body 
bound state (Efimov effect). Note, however, that in the three-body system 
the binding energy of the subsystem is not a fixed quantity but depends on 
the cm. momentum of the subsystem, i.e. the binding energy varies in the 
three-body bound state and therefore new bound states are not expected to 
appear. 

3 Reaction time scales 

One consequence of this strongly enhanced cross section can be directly noticed 
in the time scales involved. To-this end we linearize the Boltzmann equation 
given in (uh as e.g. done in Ref.H in the context of the Green function approach. 
Linearizing with respect to small variations of the deuteron distribution 5f = 
f° — f from the equilibrium one / leads to 

£r6f d (P,t) = J——Sf d (P,t), (8) 

at T hu [P,n,T\ 

where we have introduced the deuteron break-up time t^ u [P, n, T] that depends 
on the density n, temperature T, and deuteron momentum P. Using the cross 
section and the relative velocities the deuteron break-up time may be written 
as 

i- 4 t* 
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The numerical result is shown in Fig. 0, where we compare the use of the 
isolated to the in- medium break-up cross section in Eq. (@) at T — 10 MeV 
and a nuclear density of n = 0.007 fm -3 ~ no/25, no normal nuclear matter 
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Figure 3: Deuteron break-up time as a func- Figure 4: Chemical relaxation time at T ■- 
tion of the deuteron momentum. 10 MeV as a function of nuclear density. 



density. The medium dependent elementary cross section leads to shorter 
fluctuation times. As expected the difference is vanishing for larger momenta. 
A similar result is found for the relaxation time to chemical equilibrium 
as explained in the following. Linearizing the corresponding rate equations 
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with respect to small distortions 8n = n a — n from the equilibrium distribution 
n , viz. 
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In Fig. y we show the comparison of using the different cross sections under 
discussion. Again JJic difference in time scales is significant, increasing with 
increasing densityEj. 

Another well known many-body feature of bound states is the Mott effect. 
For fermions bound states should vanish beyond a certain density. This density 
is momentum dependent (Mott momentum) that is shown in Fig. pi Below 
the lines shown no pole exist in the respective few-body equations. Using 
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Figure 5: Mott momentum for deuteron Figure 6: Total number of deuterons in 
and tritons at T = 10 MeV (solid lines) and time. Due to statistics and temperature the 
T = 20 MeV (dashed lines). lines are considered preliminary. 



the present approach this is elaborated in detail for the three-body system in 
Ref.EJ. The Mott effect is an important ingredient in microscopic simulations 
of heavy ion collisions, see e.g. El, and may have consequences for the triton to 
helium-3 ratio in asymmetric nuclear matter. 



4 Extension to heavy ion collisions 

So far the investigation has been restricted to nuclear matter (homogeneous 
and symmetric). A laboratory situation closest to nuclear matter is provided 
by heavy ion collisions. As a test case we investigate the collision 129 Xe+ 119 Sn 
at 50 MeV/A lab. energy. jEhis collision has been investigated by the INDRA 
collaboration, see e.g. Ref. E2. 

The density profile suggests that the medium effects influence the deuteron 
formation in the final stage of the heavy ion collision, because of the Mott ef- 
fect. The temperature in this final stage extracted from the BUU simulation 
assuming a Fermi distribution is rather homogeneous and T ~ 4 ... 6 MeV. 
This is compatible with results we_achieved with the Quantum Molecular Dy- 
namics code provided by Aichclinll3 and other calculationsli-1 Fig. refers to 
a preliminary calculation using the .cross sections shown in Fig. that have 
been implemented in the BUU code El. The number of deuterons is increasing 
up to 10%. Also, note that the in-medium effect becomes visible at the final 
stage only. The density is small enough for deuterons to significantly survive. 

Finally, we address the question of chemical equilibrium during the heavy 
ion collision. Fig. R shows the total numbers of deuterons and nucleons as 
given in the BUU simulation (solid line) compared to the law of mass action 
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Figure 7: Number of nucleons and 

deuterons in the BUU simulation (solid 
lines) compared to the respective law of mass 
action (dashed lines). 
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Figure 8: Deuteron to nucleon ratio for the 
BUU simulation (solid line) compared to the 
law of mass action (dashed line). 



result. The upper figure ignores the Mott effect (that would lead to a wrong 
description of the experimental data) , whereas the lower figure shows the same 
with the Mott effect properly taken into account. The calculation is done 
without medium modifications for simplicity. Including the Mott effect reduces 
the number of final deuterons by about a factor of 3 and in turn increases the 
number of nucleons by about a factor of two. In the deuteron to nucleon ratio 
this ends up to about a factor of 6 shown in Fig. g for the same process. It 
seems that for the process considered ( 129 Xe+ 119 Sn at 50 MeV/A lab. energy 
central collision) the evolution is rather close to chemical equilibrium provided 
the Mott effects are properly accounted for. 



5 Conclusion and Outlook 

In summary, we have presented a systematic and consistent decoupling scheme 
to handle the many-body problem. This leads in turn to in-medium few-body 
equations that have been rigorously solved for the three-body system to treat 
the deuteron formation via three-particle collisions. The prominent changes 
induced by the medium are the self-energy shifts, the Pauli blocking, the Mott 
effect and the change of reaction rates. Within linear response we have shown 



that the change in the rates leads also to faster time scales. The deuteron 
production for the BUU simulation of the 129 Xc + 119 Sn at 50 MeV/A lab. 
energy central heavy ion collision is enhanced by 10%. 

The exact treatment of few-body systems embedded in the medium is 
however not restricted to nuclear physics. Further potential applications are 
possible, e.g., in the field of semiconductors to treat the formation of excitons 
or the trion bound state in_the dense low dimensional electron plasmas. For 
recent experiments see e.g. 113. For highly ionized dense plasmas the impulse 
approximation may fail for the three-particle break-up cross section needed to 
calculate the ionisation and recombination rates. 

Finally, the formalism is capable to also treat the four-particle problem. 
In this context the a-particle is the most interesting object. Because it is light 
and has a large binding energy per nucleon, it should play a dominant role in 
the decompressing hot nuclear matter. 
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